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Abstract 



t—( In our previous researches, a new digital watermarking scheme based 

^- on chaotic iterations has been introduced. This scheme was both stego- 

OQ secure and topologically secure. The stego-security is to face an attacker 

UN in the "watermark only attack" category, whereas the topological security 

concerns other categories of attacks. Its Lyapunov exponent is evaluated 
>— * here, to quantify the chaos generated by this scheme. 

en 

^^ Keywords: Lyapunov exponent; Information hiding; Security; Chaotic 

s ! iterations; Digital Watermarking. 

• rH 1 Introduction 

x 

J~j It currently exists only three data hiding schemes being both stego-secure and 

topologically secure. The first one is the "Natural Watermarking" with param- 
eter r) = 1 U- The two others are based on chaotic iterations. The first of them 
is a one bit watermarking scheme [71 [2] , whereas the last one allows stegano- 
graphic operations 6 . In order to enlarge the knowledge about the security 
of these processes, the Lyapunov exponent of the digital watermarking scheme 
based on chaotic iterations is evaluated here. 

This document is organized as follows. In Section [2j some basic reminders 
are given. The semiconjugacy allowing the exponent evaluation is described 
in Sect. [3| In the next one, the exponent is evaluated. This paper ends by a 
conclusion section where our contribution is summarized. 



2 Basic Reminders 

2.1 Chaotic Iterations and Watermarking Scheme 

Let us consider a system with a finite number N € IN* of cells, so that each cell 
has a boolean state. A sequence which elements belong into [1; N] is a strategy. 
Finally, the set of all strategies is denoted by [1, N] 1 *. Let S n denotes the n th 
term of a sequence S, and Vj the i" 1 component of a vector V. 

Definition 1 The set B denoting {0, 1}, let / : B N — > B N be a function and 
S £ [1, N] K . The chaotic iterations are defined by x° £ B N and 

{ T «-l jf an i ■ 

(/ (a; «-i)) if5» = [ 

In other words, at the n th iteration, only the S n — th cell is "iterated". Let 
us now recall how to define a suitable metric space where chaotic iterations are 
continuous |3]. 

Let 5 be the discrete boolean metric, S(x, y) = <=> x = y. Given a function 
/, define the function: 

F f : [1; N]xB N — ► B N 

(k,E)^(E r 6(k,j) + f(E) k .5(kJ)) 

Consider the phase space X = [1; N] w x B N , and the map defined on <Y by: 

G f (S,E) = (a{S),F f (i(S),E)), (1) 

where a : (S") nGH G [1, N] K — >■ (S n+1 )„ eW G [1, N] w and i : (5") neM € 
[1, Nj m — > S° G [1, N] are respectively the shift and the initial functions. 
Then chaotic iterations can be described by the following discreet dynamical 
system: 

{ X ° G X (2) 

J X k+1 = G f (X k ). [Z) 

To study whether this dynamical system is chaotic [S], a distance between 
X = (S, E) and Y = (S, E) € X has been introduced in 3J as follows: d(X, Y) = 
d e {E,E)+d s {S,S), where: 

N q °° I ok _ ck I 
d e (E,E) =Y6(E k ,E k ) and d s (S,S) - ^^^ ^ 



fc = l fe=l 

This distance has been introduced to satisfy the following requirements. If 
the floor value [d(X, Y)\ is equal to n, then the states E and E differ in n cells. 
In addition, its floating part is less than 10~ fc if and only if the first k terms 
of the two strategies are equal. Moreover, if the k th digit is nonzero, then the 
k th terms of the two strategies are different. With this metric, and the boolean 
vectorial negation /o, it has been proven in [3] that, 



Theorem 1 Gf is continuous and chaotic in [X, d). 

The digital watermarking scheme proposed in [71 [5] is simply the iterations 
of this dynamical system on the least significant coefficients of the considered 
media. Each property exhibited by the dynamical system will then be possessed 
too by the watermarking scheme. For further explanations, see [2]. 

2.2 The Lyapunov Exponent 

Some dynamical systems are very sensitive to small changes in their initial 
condition, which is illustrated by both the constants of sensitivity to initial 
conditions and of expansivity [3J. However, these variations can quickly take 
enormous proportions, grow exponentially, and none of these constants can illus- 
trate that. Alexander Lyapunov has examined this phenomenon and introduced 
an exponent that measures the rate at which these small variations can grow: 

Definition 2 Given / : R — > R, the Lyapunov exponent of the system com- 



1 
posed by i° £ 1 and x n+l = f(x n ) is defined by A(x ) = lim — } In /' ( : 



x 1 - 1 



Consider a dynamic system with an infinitesimal error on the initial condition 
Xq. When the Lyapunov exponent is positive, this error will increase (situation 
of chaos), whereas it will decrease if A(ajo) ^ 0. 

Example 1 The Lyapunov exponent of the logistic map [T| becomes positive 
for /i > 3, 54, but it is always smaller than 1. The tent map [9] and the doubling 
map of the circle [5] have a Lyapunov exponent equal to ln(2). 

To evaluate the Lyapunov exponent of our digital watermarking scheme, 
chaotic iterations must be described by a differentiable function on R. To do 
so, a topological semiconjugacy between the phase space X and R must be 
written. 

3 A Topological Semiconjugacy 

3.1 The Phase Space is an Interval of the Real Line 

3.1.1 Toward a Topological Semiconjugacy 

We show, by using a topological semiconjugacy, that chaotic iterations on X 
can be described as iterations on a real interval. To do so, some notations and 
terminologies must be introduced. 

Let <Sn = [1; N] W be the set of sequences belonging into [1; N] and X N = 
5n x B N . In what follows and for easy understanding, we will assume that 
N = 10. However, an equivalent formulation of the following can be easily 
obtained by replacing the base 10 by any base N. 



Definition 3 The function ip : S 10 x B 10 ->• [0, 2 10 [ is defined by: 

ip: 4=5ioxB 10 — ► [0,2 10 [ 

((S ,S\...);(£ ,...,£ 9 )) — > ¥>((£,£)) 

where (S, £) = ((S 10 , 5 1 , . . .); (-Bo, • • • , E 9 )), and ip ((S, E)) is the real number: 

9 
• whose integral part e is \2 9 ~ fe -Efc, that is, the binary digits of e are 

fc=0 
Eq El ... Eg. 



• whose decimal part s is equal to s = 0, 5° S 1 S 2 ... = J2k=i 10 _fc 5 ft - 1 . 

</3 realizes the association between a point of Afio and a real number into 
[0, 2 10 [. We must now translate the digital watermarking process G/ based 
on chaotic iterations on this real interval. To do so, two intermediate functions 
over [0, 2 10 [ denoted by e and s must be introduced: 

Definition 4 Let x e [0, 2 10 [ and: 

9 

• eo, . . • , eg the binary digits of the integral part of x: \x\ — N 2 e^. 

• (s fc )feeiN the digits of x, where the chosen decimal decomposition of x is the 

+00 

one that does not have an infinite number of 9: x = \x\ + 2_, s 10~ . 

fe=0 

e and s are thus defined as follows: 

e: [0,2 10 [ — ► B 10 

x 1 — > (e ,...,e 9 ) 

and 

*: [0,2 10 [ — ^ [0,9] M 

x 1 — > (s k )kew 

We are now able to define the function g, whose goal is to translate the 
chaotic iterations Gj on an interval of B,. 

Definition 5 g : [0, 2 10 [ — > [0, 2 10 [ is by definition such that g(x) is the real 
number of [0, 2 10 [ defined bellow: 

• its integral part has a binary decomposition equal to e , • • • , e 9 , with: 

e(x)i if % 7^ s 

e(x)i + 1 (mod 2) if i = s° 

whose decimal part is s(x) , s(x) 2 , . . . 



9 +oo 

In other words, if x = ^ 2 9 ~ k e k + ^ s k lO - * -1 , then: 

fc=0 fc=0 

9 +00 

ff(x) = ^ 2 9 - fc (efe + *(*, 8°) (mod 2)) + ^ s k+1 10- k ~ 1 . 

fc=0 fe=0 

3.1.2 Defining a Metric on [0,2 10 [ 

Numerous metrics can be defined on the set [0,2 10 [, the most usual one being 
the Euclidian distance A(x,y) — \y — x\ 2 . This Euclidian distance does not 
reproduce exactly the notion of proximity induced by our first distance d on 
X . Indeed d is richer than A. This is the reason why we have to introduce the 
following metric: 

Definition 6 Given x,y £ [0, 2 10 [, D denotes the function from [0, 2 10 [ to 
R + defined by: D(x,y) = D e (e(x),e(j/)J + D s [s(x),s(y)j, where: 

9 °° \a k — S k \ 

D e (e,e) = ^<5(e fc ,e fe ), and D s (s,s) = ^2 — ^ — • 
fc=0 fe=i 

Proposition 1 D is a distance on [0, 2 10 [. 

Proof The three axioms defining a distance must be checked. 

• D ^ 0, because everything is positive in its definition. If D(x,y) = 0, 
then D e (x,y) — 0, so the integral parts of x and y are equal (they have 
the same binary decomposition). Additionally, D s (x,y) — 0, then Vfc G 
IN*, s(x) k = s(y) . In other words, x and y have the same k— th decimal 
digit, Vfc £ IN*. And so x = y. 

• D{x,y) = D(y,x). 

• Finally, the triangular inequality is obtained due to the fact that both S 
and \x — y\ satisfy it. 

The convergence of sequences according to D is not the same than the usual 
convergence related to the Euclidian metric. For instance, if x n —^ X according 
to D, then necessarily the integral part of each x n is equal to the integral part 
of x (at least after a given threshold), and the decimal part of x n corresponds to 
the one of x "as far as required". To illustrate this fact, a comparison between 
D and the Euclidian distance is given Figure [T] These illustrations show that D 
is richer and more refined than the Euclidian distance, and thus is more precise. 



3.1.3 The Semiconjugacy 

It is now possible to define a topological semiconjugacy between X and an 
interval of R: 

Theorem 2 Chaotic iterations on the phase space X are simple iterations on 
R, which is illustrated by the semiconjugacy given bellow: 



Sio x 



? in 



G 



/o 



( [0,2 10 [,^ 



S w xB 10 ,d 



([O^f, 



D 



Proof tp has been constructed in order to be continuous and onto. 
In other words, X is approximately equal to [0, 2 N [. 



3.2 Chaotic Iterations Described as a Real Function 

It can be remarked that the function g is a piecewise linear function: it is linear 



on each interval having the form — , 

5 [10 10 

equal to 10. Let us justify these claims: 



G [[0; 2 10 x 10] and its slope is 



Proposition 2 Chaotic iterations g defined on R have derivatives of all orders 
on [0, 2 10 [, except on the 10241 points in I defined by I — / n G [0; 2 10 x 10] 



Furthermore, on each interval of the form. 



.10' 10 

g is a linear function, having a slope equal to 10: Vx ^ I,g'{x) = 10 



n n + 1 



, withne [[0;2 10 xl0] 



Proof Let I n = 



n n + 1 



, with n G [0; 2 10 x 10]. All the points of I n have 

the same integral part e and the same decimal part s°: on the set /„, functions 
e(x) and x m> s(x)° of Definition Pfl only depend on n. So all the images g(x) of 
these points x: 

• Have the same integral part, which is e, except probably the bit number 
s . In other words, this integer has approximately the same binary de- 
composition than e, the sole exception being the digit s° (this number is 
then either e + 2 10 ~ s or e — 2 10 ~ s , depending on the parity of s°, i.e., it 
is equal to e + (-l) s ° X 2 10 - s °). 

• A shift to the left has been applied to the decimal part y, losing by doing 
so the common first digit s . In other words, y has been mapped into 
10 x y-s°. 



To sum up, the action of g on the points of / is as follows: first, make a 
multiplication by 10, and second, add the same constant to each term, which is 



^ (e + (-!)•■ x*°- °) 



10 

Remark 1 Finally, chaotic iterations used in our watermarking scheme are 
elements of the large family of functions that are both chaotic and piecewise 
linear (like the tent map 0). 

We are now able to evaluate the Lyapunov exponent of our digital water- 
marking scheme based on chaotic iterations, which is now described by the 
iterations on R of the g function introduced in Definition [5j 



4 Evaluation of the Lyapunov Exponent 

Let C = ix° e [0, 2 10 [ / Vn e IN, x n £ l\, where I is the set of points in the 

real interval where g is not differentiable (as it is explained in Proposition^]). 
Then, 

Theorem 3 Vx° G C, the Lyapunov exponent of chaotic iterations having x° 
for initial condition is equal to A(x ) = ln(10). 

Proof It is reminded that g is piecewise linear, with a slop of 10 (g'(x) = 10 
where the function;; is differentiable). ThenVx € C, A(x) — lim„_j. +00 — Y^i=i m 



1 



1 



lim^+oo - Yh-i ln l 10 l = lim «^+oo -nln |10| = In 10. 
n n 



Remark 2 The set of initial conditions for which this exponent is not calculable 
is countable. This is indeed the initial conditions such that an iteration value 

Ti 

will be a number having the form — , with n 6 IN. We can reach such a real 

number only by starting iterations on a decimal number, as this latter must have 
a finite fractional part. 

Remark 3 For a system having N cells, we will find, mutatis mutandis, an 
infinite uncountable set of initial conditions x° £ [0; 2 N [ such that X(x°) — 
ln(N). 

So, it is possible to make the Lyapunov exponent of our digital watermarking 
scheme as large as possible, depending on the number of least significant coef- 
ficients of the cover media we decide to consider. Obviously, a large Lyapunov 
exponent make it impossible to achieve the well-known Original Estimated At- 
tacks @]. 



(O 



5 Conclusion and Future Works 

As a conclusion, we have available to us now a new quantitative property con- 
cerning our digital watermarking scheme based on chaotic iteration: its Lya- 
punov exponent is equal to ln(N), where N is the number of least significant 
coefficients of the cover media. This exponent allows to quantify the amplifi- 
cation of the ignorance on the exact initial condition (the media without wa- 
termark) after several iterations of the watermaking process. It illustrates the 
disorder generated by iterations of our watermarking process, reinforcing its 
chaotic nature. 

Using the scmiconjugacy described here, it will be possible in a future work 
to compare the topological behavior of chaotic iterations on X and R, and to ex- 
plore the topological security of the watermarking scheme using this new topol- 
ogy. Finally, an analogue study of the two other topologically secure schemes 
will be also conducted in order to compare these processes, being thus able to 
choose the best one according to the type of applications under consideration. 
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(a) Function x — > dist(x; 1, 234) on the interval (0; 5). 




(b) Function x — > dist(x;3) on the interval (0;5). 



Figure 1: Comparison between D (in blue) and the Euclidian distance (in green). 
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